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ABSTRACT 

An analytical model of a blister type specimen for evaluating adhes- 
ive bond strength is developed. Plate theory with shear deformation Is 
used to model the deformation of the plate^and elastic deformation of the 
adhesive layer Is taken Into account. It Is shown that the Inclusion of 
the elastic deformation of the adhesive layer can have a significant In- 
fluence in the energy balance calculations of fracture mechanics. 

INTRODUCTION 


The increasing use of adhesives In the development of high strength 
composite materials and in the joining of structural components has brought 
about the need for standard methods of testing the strength properties of 
adhesives. One such test proposed for use under a variety of environmen- 
tal conditions involves a blister specimen subjected to pressure loading 
[1]. [2]. This test employs a circular plate, usually of uniform thick- 
ness, bonded through a concentric annular adhesive layer to a flat sup- 
porting surface. During testing, the plate Is loaded by means of fluid 
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pressure applied to the unbonded concentric center portion or the plate. 
Figure 1 shows a sketch of a circular plate blister specimen of radius b 
bonded over the annular area a<x<b to the rigid supporting surface. 

After a test to Initial debonding or to failure, the specific adhesive 
fracture energy y. of the bond Is calculated by means of an energy balance* 

a 

which relates y. to the pressure loading at Initiation of debonding and 

O 

the elastic and geometric properties of the specimen. 

The curves given In [2] for determining y. are calculated by either 

O 

plate theory or the finite element method. In these calculations It Is 
assumed that the adhesive layer Is so thin that the plate may be consid- 
ered to be bonded rigidly to the support neglecting the compliance of the 
adhesive layer. Such an assumption is reasonable when the elastic modu- 
lus of the plate has the same order of magnitude as that of the adhesive 
and the thickness of the adhesive layer Is small in comparison with the 
plate thickness. However, when the circular plate Is metal, the adhesive 
can have a much lower elastic modulus. In such a case the compliance of 
the adhesive layer may no longer be neglected In the energy balance cal- 
culations. 

The purpose of this paper Is to Investigate the effect of the adhes- 
ive layer compliance on the fracture mechanics of thin plate blister 
specimens. Since linear plate theory Is employed, the results are rn~ 
stricted to small values of h/a. Also, the maximum plate deflection must 
be somewhat less than the plate thickness. Because of the similarities 
of analysis of circular plate bending and cylindrical bending under plane 
strain conditions, both cases are treated In this paper. 
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ANALYSIS 


The present analysis treats the plate using equations of a plate theory 
which Includes deformations of bending* In-plane extension, and transverse 
shear. For the one dimensional problems of cylindrical bending (plane 
strain) or of axlsymmetrlc deformation of a plate the displacement field 
Is described completely by means of the variables u, w, and 6 . The stres- 
ses In the plate are related to the stress resultants and stress couples 
I'll* N 2 * Q, H|» and M 2 through the usual formulas of plate theory, where 
the In-plane normal stress varies linearly across the thickness, the 
transverse shearing stress varies parabollcally, and the transverse nor- 
mal st'i'ess Is considered to be small. The Independent variable x Is 
shown In Figure 1. Plate theory equilibrium equations are 


dN^/dx » T + c.|(N 2 - N^)/x 

( 1 ) 

dQ/dx ■ -q - c^ Q/x 

( 2 ) 

dM^/dx » m + Q + c^(M 2 - M^)/x 

(3) 


where c.j Is equal to 0 for the plane strain case and equal to 1 for the 
axlsymmetrlc case. The distributed moment m Is given by 

m * h T (A) 

where h us equal to h /2 If coupling between bending and extension Is con- 
sidered and Is 0 If this effect Is neglected. The strain displacement 
equations are 

62 • c^u/x , du/dx = e-j (5) 
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dw/dx ■ e - Y 


( 6 ) 


• c^e/x , d3/dx • Ky (7) 

These equations are coupled through the stress-strain relations 

■ K(e^ + VC2) (8) 

N 2 » K(c 2 + ve^) (9) 

Y - + m/6)/C (10) 

= D(ic.| + VK 2 ) (^1) 

M 2 » D(ic2 + VK-|) 02) 

where the plate rigidities K, C, and 0 are given by 

K=Eh/(l-v^) (13) 

C « 56h/6 « 5(1 - v)K/12 (14) 

0 « Kh^/12 (15) 


and the constant C 2 Is to be set equal to 1 If transverse shear deforma- 
tion Is Included or set equal to 0 If It Is Ignored. 

The thin adhesive layer has an elastic modulus smaller than 
that of the plate; therefore. It Is reasonable to Ignore the In-plane 
stress components In the adhesive layer when computing equilibrium of the 
plate and adhesive layer together. Hence, the adhesive layer Is treated 
as a distributed spring or elastic foundation which transmits normal and 
shear stresses between the plate and the rigid support. Accordingly, 
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the surface loads on the bonded portion of the plate are related to the 
plate displacements by means of the expressions 

q • 0 • kjjw ( 16 ) 

for the normal component and 

T « kj(u ♦ h0) (17) 

for the shear components. The corresponding surface loads on the unbonded 
portion of the plate are 

q « - p (18) 

and 

T • 0 ( 19 ) 

respectively. 

The elastic foundation rigidity constants are 

St • Vo 

for the normal stress and 

S ■ V'o 

for the shear stress. The elastic constant Is equal to If the nor- 

mal stress Is assumed to be uniaxial* Is equal to E^/(l - v^) If one In- 
plane stress and one In-plane strain vanish* and Is equal to (1 - ^q)6q/ 

(1 + Vjj)(l - 2 v^) If both In-plane strains are assumed to be zero. 

If coupling of bending and extension Is neglected (h«0) and trans- 
verse shear deformation Is Ignored (03 ” 0 )* the governing equations of 
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the plane strain problem reduce to 


w'"' » P/D (22) 

for the unbonded portion of the plate and to 

(23) 

In the bonded portion. Here, primes denote differentiation with respect 
to X. The length which characterizes the stiffness ratio of the plate 
to the elastic foundation. Is given by 

I - (D/k,,)^ (24) 

The general solutions of (22) and (25) are 

w ■ A.J + A 2 X 

+ AjX^ + A^x^ + px^/24D (25) 

and 

w ■ Age®’^ cosax + Age®’* sinew 

Aye"®’* cosox + Age"®’* slnox (26) 

respectively, where the quantities A^ are arbitrary constants and 
a ■ ^/t . 

The constants A^ are found from boundary conditions 
3*0 • Q»0 , atx*0 
M^*0 , Q*0 , atx*b 
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dnd fliatchlng conditions of w» B» and Q at x ■ a. If (b - 3, 

as would ba the case for many practical blister test specimens* the con- 
stants Ag and Ag become negligible and the boundary conditions at x*b 
have negligible Influence on the solution. 

Once the solution for w(x) has been found* one can evalua*" the max- 
imum normal stress In the adhesive layer using 

®max * *^*^^*^ 

Following through the calculations one obtains the expression 
where 

- (1 + 3v7 £/a + 3£W)/(1 + v7-t/a) 

For the axlsymmetrlc problem the system of equations reduces to 

x"^{x[x"^xw' P/D (29) 

for 0 < X < a and to 

x“^{x[x'^(xw’ )•]•}' - -w/f^ (30) 

for a < x < b when h * 0 and c^ • 0. 

The general solutions of (29) and (30) are [3] 

O 

w ■ B^x logx + B 2 logx 

+ Bgx^ + B4 ♦ pxVsAD (31) 
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and 


w ■ Bg ber(x/<C) *■ Bg ba1(x/i) 

♦ By ker(x/i) ♦ Bg ke1(x/^) (32) 

respectively* where her* be1» ker* and kel are zero order Kelvin functions. 

The boundary conditions that B and Q vanish at x • 0 require that 
B.| ■ 0 and B 2 ■ 0. The other constants are found from the boundary con- 
ditions 


Hj»0 * Q«0 atx»b 

and matching conditions of w» B* M^» and Q at x ■ a. If (b - i)it > 3, 
the constants Bg and Bg may be set equal to zero, and the plate may be 
treated as Infinite; hence, the length b will not appear In the solution. 
Also, If a/£ » 1* some simplification Is achieved by using one term 
asymptotic formulas for the Kelvin factions. It Is expected that these 
conditions will be satisfied for many practical blister test specimens. 

Calculations of the maximum normal stress In the adhesive layer show 
that Its value Is given by 

®max^P • (a^/8t^)ti)g (33) 

where 

- [fj - 4(C/a)f, ♦ - 2(^/a)f4] 

f^ « ker^(a/£) + ke1^(a/^) 
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f 2 • ker(a/t)ke1'(a/t) - ke1(aM)ker'(aM) 


fj • ker(a/£)ker'(a/4) ♦ ke1(aM)ke1'(a/£) 

• ker'2(a/^) + kei'^(aM) 

(b - a)M j 3 

For cases where coupling of extension and bending effects Is Included 
and transverse shear deformation Is considered* numerical solutions to the 
governing system of differential equations and boundary conditions may be 
obtained by means of any one of the numerical Integration techniques which 
have been applied successfully to the analysis of shells of revolution. 
Results reported later In this paper were obtained by using a multl- 
se^nent numerical Integration technique [4] to calculate the plate dis- 
placements. stress results, stress couples, and elastic foundation reac- 
tions at discrete x values. 

CALCULATION OF 3U/3A 

Bennett. Devries, and Williams [1] have proposed the adhesive frac- 
ture criteria 

3U/3A » Yj (34) 

for a blister test specimen under pressure loading. Here U Is the total 
strain energy In the elastic system and A Is the debonded area, which Is 
equal to the length a for a plane specimen of unit width and equal to 
ira^ for the circular specimen. The parameter y. Is the specific adhes- 
ive fracture energy, representing the energy required to debond a unit 


area. 
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For a linearly elastic system the strain energy U may be computed 
from values of the pressure loading and the displacements of the pressur* 
Ized portion of the plate according to 


2U • 


f 

Jo 


pw dx 


(35) 


for the plane problem and 

2U • 2ir [ pM X dx 
^0 


(36) 


for the axisymmetric problem. The displacement w must be found by solving 
a boundary value problem. In [1]» classical plate theory for thin speci- 
mens and a finite element method for thicker specimens are used to calcu- 
late 3U/3A for a plate bonded rigidly to a rigid support; I.e. the com- 
pliance of the adhesive layer Is neglected. Using deflection expressions 
of classical plate theory the values of 3U/3A turn out to be 


3U/3A - pV/18D 


(37) 


for the plane specimen and [1] 
3U/3A • pV/128D 


(38) 


for the circular plate. In order to apply the fracture criterion (34) to 
specimens with compliant adhesive It Is necessary to calculate 3U/3A for 
the plate on elastic foundation model. 

Although the evaluation of Integrals (35) and (36) for U and the diff- 
erentiation with respect to A could In principle be carried out numerically 
for the plate on an elastic adhesive layer, a much simpler method for com- 
puting 3U/3A for this case does exist. Its derivation Is given In the 
Appendix. 
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Closed form expressions for dU/3A may be obtained from (28) and (33) 
together with (A-8) of the Appendix (t ■ 0). The results are 

3U/3A • (w, ♦ 3t^/a^)^pV/18D (39) 

for the plane strain case and 

3U/3A • (wg ♦ 8£^/a^)^oV/128D (40) 

for the axl symmetric case. Here and u>2 are functions of t/a defined 
previously In (28) and (33). 

RESULTS AND DISCUSSION 

Results of calculations of using (28) and (33) are given In 

Figure 2. The effect of the adhesive layer compliance on the calculated 
value of Is Indicated by the dependence on the dimensionless para- 
meter t/i. For small values of C/a the values of reduce to 

• (a«)*/n (41) 

where n • 8 for the axl symmetric case and n ■ 3 for the plane strain case. 
In terms of given dimensions and elastic constants of the plate and ad- 
hesive layer the ratio C/a Is 

C/a • [E/(l - v2)BQ]^[h\/12a^]^ (42) 

where the first factor Indicates the Influence of material properties and 
the second Indicates the Influence of geometry. From (41) and (42) one 
can see that becomes singular as h^-^D, the order of the singular- 

ity being given by 
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of the differential equations were obtained for aluminum plate-epoxy ad- 
hesive test specimens. Figure 4 shows typical curves for the stress dis- 
tribution In the adhesive layer. Results of other cases are summarized 
In Table 1. In Table 1 the values enclosed In parantheses are based on 
the closed form solutions. Numerical analysis for cases with no coupling 
of bending and extension and without plate shear deformation show that 
the numerical analysis method produces results within one half percent of 
the closed form solutions. It appears that for specimens In the range 
Included In Table 1 the cloi^d form solutions for the ratio are 

about 10% too high. HOk.«». * , the Inclusion of these second order effects 
appears to make negligible change In the calculated values of dU/dA, even 
for a ratio h/a as high as 0.25. This does not mean that the results of 
the closed form solutions are to be considered accurate for thick plates. 
For thick plates the hypothesis of plane sections, assinned in all the 
theories employed In this paper, may not be satisfied with sufficient 
accuracy. In fact, finite element solutions of [1] Indicates that con- 
siderable deviations from this hypothesis exist for plates with rigid 
adhesive. Calculations show that the ratio lies in the range 

0.20 - 0.25 for all specimens reported in Table 1. 

Consider now how the material properties of the plate and adhesive 
layer enter into the parameter £/a. From (42) it Is seen that the material 
factor Is [E/(l - where the elastic constant for the adhesive 

layer depends on the Young's modulus, the Poisson's ratio, and an assump- 
tion regarding the in-plane stress or strain components within the layer. 
Of the possible in-plane stress or strain conditions to be assumed, it 
would seem that the assumption of no In-plane strain would be most 
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appropriate where a low modulus adhesive Is used together with a high mod- 
ulus plate and support. In such a case, the normal strain In the adhesive 
layer Is much larger than the In-plane components. Of course, this assump- 
tion Is not valid right at the edge of the bond region where one of the 
stress conditions might be more appropriate; however, at distances on the 
order of h^ away from the edge. It is expected that In-plane constraint 
of the adhesive layer Is achieved within reasonable accuracy. 

For the case of In-plane constraint of the adhesive, the material 
factor of the ^/a ratio becomes 

[E/O - » Ic(E/Eq)^ 

where 

k » [(1 + Vq) 0 - 2v^)/(l - v2)(l - 

The factor k Is not very sensitive to Poisson's ratio v of the plate; 
however It becomes quite sensitive to Poisson's ratio of the adhesive 
for values of close to 0.5 as can be seen In Table 2. 

Since the material factor appraoches zero for an Incompressible ad- 
hesive, t/a also approaches zero. For the case t/a « 0, the maximum 
normal stress In the adhesive layer becomes singular with the dis- 
tributed reaction of the elastic foundation model of the adhesive layer 
becoming a concentrated line load at the edge. However, since there Is 
no constraint of the shear deformation of the adhesive layer, this layer 
continues to behave as a distributed shear spring with finite values 


14 



CONCLUSIONS 


A structural mechanics approach has been used to Investigate the 1n> 
fluence of the elasticity of the adhesive layer on the fracture mechanics 
parameters of a blister test specimen. It has been demonstrated that In 
some practical cases the effect of adhesive layer elasticity cannot be 
neglected. A parameter which serves to Indicate the Importance of ad- 
hesive layer compliance In practical calculations has been successfully 
Identified and Is presented In terms of the product of a material factor 
and a geometry factor. 

While the efforts of this paper have concentrated on the blister 
test specimen under uniform pressure loading, the methods of analysis 
presented can also be applied In the case of a blister specimen subjected 
to displacement loading. 

Some shortcomings of the present structural mechanics approach are 
Its Inability to determine an upper limit of the ratio h/a for which the 
results are valid and Its Inability to determine details of the stress 
distribution at points near the edge of the bond region. For this reason 
a theory of elasticity study of a blister test specimen with a compliant 
adhesive layer would be desirable. With these limitations in view, it Is 
hoped that the results presented here will prove useful In the process of 
standardization of blister specimens for adhesive fracture testing. 
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NOMENCLATURE 


a 

A 

b 

®o 

C 

D 

E 


"o 

ti 

Si 

‘‘s 

K 

I 


M,.M2 

N^.Ng 

P 

q 

Q 

u 


Inner dimension (see Figure 1) 
debonded area 

outer dimension (see Figure 1) 

elastic constant for adhesive layer 

constants equal to 0 or 1 

transverse shear stiffness of plate 

bending stiffness of plate 

Young's modulus of plate material 

Young's modulus of adhesive layer 

modulus of rigidity of adhesive layer 

plate thickness 

thickness of adhesive layer 

half thickness of plate or zero 

adhesive layer transverse normal stiffness 

adhesive layer shear stiffness 

extenslonal stiffness of plate 

characteristic length 

moment of surface loading on plate 

plate bending moment components 

plate membrane stress resultants 

applied normal pressure 

normal surface loading on plate 

transverse shear resultant In plate 

In-plane displacement at middle surface of plate 
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NOMENCLATURE (continued) 


U total strain energy 

w transverse displacement of plate 

X coordinate (radius In axlsymnetrlc problems) 

6 rotation of normal In plate 

Y transverse shear strain In plate 
specific adhesive fracture energy 

e^,e2 membrane strain components In plate 
k<I»ic 2 curvature change components In plate 

V Poisson's ratio of plate material 
Poisson's ratio of adhesive layer 

T shearing stress In adhesive layer 

a normal stress In adhesive layer 

»o>2 function of t/& 
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APPENDIX 


Consider two configurations of the blister specimen subjected to the 
same value of applied pressure p. In configuration I with total strain 
energy Uj the adhesive layer acts as a continuous spring over the coordin- 
ate range a<x<b. In configuration II with total strain energy Ujj. the 
adhesive layer Is debonded over the range a<x<a * L& and acts as a con- 
tinuous spring over the range a‘t*Aa<x<b. Let AA be new area of debond- 
ing In passing from configuration I to configuration II. The derivative 
9U/3A Is then given by 


3U/3A - 11m (Uxt - U,)/AA (A-1) 

AA^O “ * 

Let [P be the elastic strain energy stored In the circular plate plui the 
strain energy stored In the continuous spring over the range ai-Aa<x<b. 
The total elastic strain energy Is then given by 

U - + (a^/2kj, + T^/2k^)AA (A-2) 

where the last term represents the elastic strain energy of the portion of 
the adhesive layer In the range a<x<a + Aa. Noting that (o»t) have 
known values configuration I and values (-p,0) In config- 

uratlon II one can write 


u„ . u, 


The strain energy U° 


«ii - “i 

Is a fiRictlon of the loads acting on the structure 
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consisting of the entire circular plate and the portion of the adhesive 
layer In the range a * Aa<x<b. Thus. 

• U'(P.F..F,) 

where and are the normal and tangential conq)onents, respectively, 
of the force transmitted to the structure at x ■ a e. These force com- 
ponents are given by 

F|(j ■ -a AA 

Fj - -T M 

Since they are carried to the structure through the portion of the adhesive 
layer which debonds In passing from configuration I to configuration II. 
Since the pressure p remains constant during debonding, the change In 
may be expressed as 

Uf, - uj • (M(“/i(F„)4F„ ♦ OU°/8Fj)4Fj (A-4) 

where the partial derivatives may be evaluated In configuration I. The 
force Increment components In passing from configuration I to configura- 
tion II are 

■ <»«x t F*** 

(A-5) 

• ■'.ax F* 

Now combining (A-1) through (A-5) and taking the limit AA**0 results In 
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W8A ■ (SU“/»Fh)(<>„x * •*> * 

♦ ■ ’«x^/2k, (A-6) 

According to Castigll arm's second theorem* the partial derivatives 3U^/9F|^ 
and are equal to the plate displacement components at the point of 

application of F|| and F^. In turn* these displacement components may be 
expressed In terms of and through the adhesive layer compliance 
relationships. Thus, 

- (A-7) 

* Ws 

Finally, substituting (A-7) Into (A-6) produces the required expression 
for 3U/3A. 

W3A . (0^ ♦ p)*/2k^ ♦ (A-8) 

It might be noted that the derivation of (A-8) did not make use of plate 
theory approximations; therefore. It is also valid when the circular 
plate Is modeled as an elastic continuum Instead of a structural plate. 
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TABLE 1. RESULTS OF NUMERICAL ANALYSIS 
(E/Ep - 22.4, V • 0.333, ■ 0.35) 


Code 

h/a 

V* 



^max/ P 

t/a 

PA.O 

.125 

.010 

1.193 (1.206) 

71.0 (78.8) 

15.6 

.071 

Pt.l 

.125 

.010 

1.222 (1.206) 

72.0 (78.8) 

16.8 

.071 

Ax.O 

.125 

.010 

1.248 (1.271) 

27.4 (30.4) 

5.90 

.071 

Ax.1 

.125 

.010 

1.290 (1.271) 

28.2 (30.4) 

6.33 

.071 

Pt.l 

.0625 

.004 

1.081 (1.097) 

282. (324.) 

76.3 

.034 

Pt.l 

.0625 

.010 

1.126 (1.123) 

192. (209.) 

40.7 

.042 

Ax.l 

.0625 

.004 

1.109 (1.123) 

108. (124.) 

28.7 

.034 

Ax.1 

.0625 

.010 

1.162 (1.156) 

74.0 (80.2) 

15.3 

.042 

Pl.l 

.25 

.010 

1.409 (1.368) 

28.8 (31.0) 

6.80 

.119 

Ax.1 

.25 

.010 

1.550 (1.482) 

11.5 (12.0) 

2.55 

.119 


Code: 91 . • Plane Strain Case 

Ax. * Axlsymmetrlc Case 

0 ■ Shear Deformation of Plate Ignored 

1 > Shear 0efom.nt1on of Plate Included 
n ■ 3 for plane strain case 

n ■ 8 for axlsynmetrlc case 
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TABLE 2. VALUES OF k. (v - 1/3) 


''o 

.25 

.3 

.35 

.4 

.45 

k 

.98 

.96 

.91 

.85 

.74 



.48 

.49 

.499 

.4999 

.5 

k 

.60 

.51 

.29 

.16 

.00 
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Fl 9 ure 1. Specimen for Blister Test of Adhesive Bond. 

Figure 2. vs. A/a. 
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and Extension only* Transverse Shear Deformation Neglected; 
Curve (c)» Bending* Extension* and Shear. 
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(a) Axlsynmetrlc case (n » 8) 

(b) Plane strain case (n » 3) 


Figure 2. 




V^DOU/dA) /a*p 



Figure 3. *^n^b(au/3A)/a*p vs. 4/a. 

(a) Axl symmetric case (n > 8) 

(b) Plane strain case (n • 3) 



T/p 



Figure 4. 


o/p and t/p vs. x/a by numerical Integration method for axl- 
symmetric case. Curve (a). Bending only; Curve (b). Bending 
and Extension on1y» Transverse Shear Deformation Neglected; 
Curve (c). Bending* Extension* and Shear. 
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